Interface-anisotropy induced asymmetry of intermixing in bilayers 

P. Sule, M. Menyhard, L. Kotis, J. Labar, W. F. Egelhoff Jr.* 
Research Institute for Technical Physics and Material Science, 
Konkoly Thege u. 29-33, Budapest, Hungary, sule@mfa.kfki.hu,www.mfa.kjki.hu/^sule, 
* National Institute of Standards & Technology, Gaithersburg, Maryland 20899 

(February 2, 2008) 

The ion-sputtering induced intermixing is studied by molecular dynamics (MD) simulations and 
by Auger electron spectroscopy depth profiling (AES-DP) analysis in Pt/Ti/Si substrate (Pt/Ti) 
and Ta/Ti/Pt/Si substrate (Ti/Pt) multilayers. Experimental evidence is found for the asymmetry 
of intermixing in Pt/Ti and in Ti/Pt. An unexpected enhancement of the injection of the heavy 
Pt atoms into the Ti substrate is observed both by AES-DP and by MD simulations. In Ti/Pt we 
get a much weaker interdiffusion than in Pt/Ti. The asymmetry is explained by the backscattering 
of hyperthermal particles at the mass-anisotropic interface and which is reproduced by computer 
atomistic simulations. The AES-DP measurements support our earlier predictions (P. Siile, M. 
Menyhard, Phys. Rev., B71, 113413 (2005)) obtained for mass-anisotropic bilayers. 
PACS numbers: 61.80. Jh, 68.35.-p, 68.35.Fx, 66.30.-h, 68.55.-a 



There has been intense interest in the past few years 
in maintaining well-defined conditions for the fabrication 
of ultra-thin film devices [1-4]. Controllable fabrication 
of interfaces is still a great challange for nanoscience and 
nanotechnology [1,3,4]. It is necessary to understand the 
driving force of intermixing in order to avoid or suppress 
interdiffusion during thin film growth. Unfortunately our 
knowledge on the mechanism of interfacial mixing is far 
from being complete [5-7] . 

Despite the apparent technological importance of thin 
film multilayers, the research of the atomic transport phe- 
nomena at interfaces, however, is less intense. Auger 
electron spectroscopy (AES) depth profiling analysis has 
been used for the study of the sharpness and broad- 
ening of interfaces during ion-sputtering [8] in combi- 
nation with transmission electron microscopy (TEM) 
[9,10]. Computer simulations have revealed that mass- 
anisotropy of the bilayer governs interdiffusion at the in- 
terface and greatly influences surface morphology devel- 
opment during ion-sputtering [11,12]. 

In this Letter we present results which can be taken as 
the first direct experimental evidence of mass effect on 
interfacial mixing. We reproduce and explain the exper- 
imentally found asymmetry of intermixing by MD simu- 
lations. 

Classical constant volume molecular dynamics simula- 
tions were used to simulate the ion-solid interaction using 
the PARC AS code [13]. The computer animations can be 
seen in our web page [14]. Further details are given in 
[13] and details specific to the current system in recent 
communications [11,15]. We irradiate the bilayers Pt/Ti 
and Ti/Pt (9 monolayers, ML film/substrate) with 0.5 
keV Ar"*" ions repeatedly with a time interval of 10-20 
ps between each of the ion-impacts at 300 K which we 
find sufficiently long time for the termination of interdif- 
fusion, such as sputtering induced intermixing (ion-beam 
mixing) [15]. The initial velocity direction of the impact- 
ing atom was 10 degrees with respect to the surface of the 
crystal (grazing angle of incidence) to avoid channeling 



directions and to simulate the conditions applied during 
ion-sputtering. We randomly varied the impact position 
and the azimuth angle 4>. In order to approach the real 
sputtering limit a large number of ion irradiation are em- 
ployed using automatized simulations conducted subse- 
quently together with analyzing the history files (movie 
files) in each irradiation steps. In this article we present 
results up to 100 ion irradiation which we find suitable 
for comparing with low to medium fluence experiments. 
100 ions are randomly distributed over a 50 x 50 area 
which corresponds to < 10^^ ion/cm^ ion fluence and the 
removal of few MLs. The size of the simulation cell is 
110 X 110 X 90 A including 57000 atoms. In order to 
reach the most efficient ion energy deposition at the in- 
terface, we also initialize recoils placing the ion above the 
interface by 10 A and below the free surface in the 9 ML 
thick film. In this way we can concentrate directly on 
the intermixing phenomenon avoiding many other pro- 
cesses occur at the surface (surface roughening, sputter 
erosion, ion-induced surface diffusion, cluster ejection, 
etc.) which weaken energy deposition at the interface. 
Further simplification is that chanelling recoils are left to 
leave the cell and in the next step these energetic and 
sputtered particles are deleted. 

The experimental setup is the following: According to 
the XTEM results the thickness of the layers in samples 
are: Pt 160 A/Ti 130 A/Si substrate, and Ta 260 A (cap 
layer to prevent oxidation of Ti)/Ti 150 A/Pt 140 A/Si 
substrate. The interfaces are sharp and straight. For 
the sake fo simplicity we consider our multilayer samples 
as bilayers and we study the atomic transport processes 
at the Ti/Pt and Pt/Ti interfaces. We see no asymme- 
try in the interface widths of Ti/Pt and Pt/Ti bilayers 
(which includes roughness and mixing) using XTEM in 
Ta/Ti/Pt/Si (Ti/Pt) and in Pt/Ti/Si (Pt/Ti) in the as 
grown samples. Both samples were AES depth profiled 
by applying sputtering conditions as following: ion en- 
ergy 500 — 1500 eV, projectile Ar+, grazing angle of inci- 
dence 5° — 10° degree (with respect to the surface), sam- 
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pie was rotated during ion bombardment. The atomic 
concentrations of Pt, Ti, Ta and Si were calculated by 
the relative sensitivity method taking the pure material's 
values from the spectra. The oxygen atomic concentra- 
tion was calculated by normalizing the measured oxygen 
Auger peak-to-peak amplitude to Ti02 [16]. The depth 
scale was determined by assuming that the sputtering 
yield (Yi) in the mixed layer is the weighted sum of the 
elemental sputtering yields (X^i^i^i)- Two typical AES 
depth profiles recorded on samples of Pt/Ti and Ti/Pt 
at 500 eV ion energy, are shown, respectively in Fig. 1. 
It is clear without any detailed evaluation that the Pt/Ti 
and Ti/Pt transitions are very different. While the Ti/Pt 
transition is a "normal" one, in case of Pt/Ti an unusual 
deep penetration of Pt to the Ti phase is observed. The 
interface widths between pure materials can be well char- 
acterized by the distance of points on the depth profile 
showing concentrations of ^ 84 % and ^ 16 % (step func- 
tion true profile) if the transitions can be approximated 
by an error function. It can be supposed that the intrin- 
sic surface roughness, Uroughi is independent of ion bom- 
bardment induced mixing and thus the measured width, 

CTmeas, Can bc Calculated as, a = (rl-^eas - ^lough- ^^'^ 
a is the width as a result of the ion bombardment in- 
duced processes (mixing) and the finite inelastic mean 
free path of the signal electrons. The metal/substrate 
system interface was found to be sharp, thus correction 
is not necessary, while roughness of the other interfaces 
are about 15 A. Though the Pt/T i transition is far from 
being an error function type one, we still used the defi- 
nition for comparative purpose. 

In Fig. 1 we also recognize the presence of oxygen in 
the Ti layer in the Ti/Pt and metal/substrate interfaces. 
The later is the consequence of using Si substrate with 
native oxide. Since the bulk level of oxygen slightly cor- 
related with the ion current intensity (the larger the ion 
current intensity the lower the oxygen AES signal) part 
of the oxygen is contamination occurring during the AES 
depth profiling process. The concentrations of the oxy- 
gen in the Ti/Pt interface neither correlate with the bulk 
level nor with the ion current intensity. Consequently we 
might conclude that the origin of the oxygen on the inter- 
face is due to the sample preparation process. Moreover 
the interface broadening does not correlate with the con- 
centration of the oxygen in the Ti/Pt interface. Thus we 
conclude that the atomic transport is not affected by the 
presence of the oxygen. 

The MD simulations of ion-sputtering and AES-DP re- 
veal rather similar values for broadening at the interface 
as well as reproduce the asymmetry of intermixing shown 
in Figs la and lb. In Fig 2 the evolution of the sum of 
the square atomic displacements of all intermixing atoms 
(i?2) ^ Y.f[Ti{t) - ri{t = 0)]2 (ri(t) is the position vec- 
tor of atom 'i' at time t), can be followed as a function of 
the simulated ion fluence, and calculated as it has been 
given in ref. [11]. A result for the Cu/Co system is also 
given to compare the magnitudes of intermixing in vari- 
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FIG. 1. The concentration depth profile as a function of the 
removed layer thickness (A) obtained by AES depth profihng 
analysis using ion-sputtering at 500 eV ion energy in Pt/Ti/Si 
substrate (upper Fig: la) and Ta/Ti/Pt/Si substrate (lower 
Fig: lb). 



ous bilayers and to demonstrate the performance of the 
MD simulations. The magnitude of ion-beam intermix- 
ing can also be characterized by the measurable quan- 
tity mixing efficiency ^ [11,18] which is given as follows: 
^ = wp^^ where $, and Fu, are the ion fluence (ion/ A ) 
and the deposited nuclear energy at the interface (eV/A), 
respectively. 

We get mixing efficiency value of ^ w 10 A /eV which 
is comparable with our experimentally measured value 
of ~ 3 — 6 A /eV for Cu/Co. The asymmetry of mixing 
can clearly be seen when (i?^)(<l>) and the depth profiles 
are compared in Ti/Pt and in Pt/Ti in Figs la, lb and 2. 
The computer animations of the simulations also reveal 
the stronger interdiffusion in Pt/Ti [14]. Moreover, in 
the inset in Fig 2 the higher mobility (oc (i?^)) of the 
energetic Pt atoms can be seen in Pt/Ti. 

A relatively weak intermixing is found in Ti/Pt (the 
corrected ion-induced AES-DP interface width a IsA 
at 500 eV ion energy) and a much stronger interdiffusion 
occurs m the Pt/Ti bilayer (cr w 70 A). MD simulations 
provide 4 (- lOA) and ~ 8 (-■ 20A) ML thick interface 
after 100 ion impacts, respectively. As mentioned before, 
in our immiscible reference system, in Cu/Co we find also 
a weak mixing. 

Asymmetric AES depth profiles (when the broaden- 
ing of A/B interface is different from that of B/A) have 
already been observed [8,19]. The asymmetric behav- 
ior was explained by the large relative sputtering yield 
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FIG. 2. The simulated {R^) (A^) in Pt/Ti, Ti/Pt and in 
Cu/Co as a function of the ion-fiuencc (ions/cm^) obtained 
during the ion-sputtering of these bilayers at 500 cV ion en- 
ergy. The dotted hne (iso Pt/Ti) denote the results obtained 
for the artificial mass-isotropic Pt/Ti bilayer. Inset: The sim- 
ulated {R^) for Pt in Pt/Ti and in Ti/Pt at 500 eV ion energy. 




FIG. 3. The schematic view of the hyperthermic atomic 
backscattering mechanism in anisotropic bilayers. The closed 
light regions represent the cascade region occurs at each ion 
impact. Note that the lower lobe of the cascade is larger 
in Pt/Ti than the upper lobe in Ti/Pt due to the downward 
scattering and the confinement of the light particles indicated 
by the curved arrows. The arrows with dashed lines show 
the atomic transport of Pt. The path of the inpinging ions 
from the vacuum towards the free surface of the film/substrate 
system is also shown by arrows. 



of the elements (preferential sputtering) causing differ- 
ent atomic mixing [8] and/or enhanced bombardment 
induced interface roughening. The asymmetry of inter- 
mixing is also known during thin film growth [6,20]. In 
the present experiment the relative sputtering yield of 
Ypt/Yn ~ 0.7 - 0.9 at 500-1500 cV therefore the mech- 
anism is different. This kind of an ion-bombardment in- 
duced asymmetric mixing in metallic diffusion couples 
has not reported yet at best of our knowledge. 

In recent papers we explained single-ion impact in- 
duced intermixing as an interdiffusion process governed 
by the mass-anisotropy parameter (mass ratio) in these 
bilayers [11,12]. Following this way of reasoning Ti/Pt 
(Pt/Ti) is a mass-anisotropic, while Cu/Co is a mass- 
isotropic bilayer. It has already been shown in refs. 
[11,15] that the backscattering of the hyperthermal par- 
ticles (BHP) at the heavy interface leads to the increase 
in the energy density of the coUisional cascade. We have 
found that the jumping rate of atoms through the inter- 
face is seriously affected by the mass-anisotropy of the 
interface when energetic atoms (hyperthermal particles) 
are present and which leads to the preferential intermix- 
ing of Pt to Ti [11,15]. In this article we would like to 
explain the mixing asymmetry seen in the case of Pt/Ti 
and Ti/Pt using the extension of this model. 

The schematic view of the hyperthermic backscatter- 
ing mechanism is shown in Fig 3. According to the 
BHP model, in Pt/Ti, the energetic Ti light particles are 
backscattcrcd downwards and confined below the heavy 
interface in the bulk which results in increasing energy 
density below the interface. The increased lifetime and 
energy density of the collisional cascade or heat spike re- 
gion below the interface in Pt/Ti leads to an enhanced 
injection of Pt atoms into the bulk and to a rather deep 
depth distribution of Pt shown in Fig la. We find a 
18 - 22 A (6 - 8 ML) deep penetration of Pt into the 



Ti phase (and 36 - 44 A interface width) at $ « 10^^ 
ion/cm^ fluence by MD simulations. Although the ex- 
perimental (J « 70A value is still far from our simulated 
value, however, we expect that the penetration depth 
(intermixing length) of Pt could be increased by further 
ion-impacts. Nevertheless, MD simulations reproduce 
the mixing asymmetry and we are able to explain the 
essential features of the phenomenon. In Ti/Pt the Ti 
atoms are backscattered in the overlayer film and ejected 
towards the surface. This leads to the reduction of the 
energy density in the collisional cascade and hence to a 
weaker intermixing when compared with the Pt/Ti sys- 
tem. The BHP is the most effective when the projected 
range of the deposited energy is positioned right at the 
interface. 

In order to clarify the details of intermixing, simula- 
tions have been run with mass ratio 5 is artificially set 
to ^ « 1 (mass-isotropic), and we find that (i?^) drops 
below the corresponding curve of Ti/Pt in Pt/Ti (see 
Fig 2, dotted line, and animation [14]). We find it as 
a direct evidence of mass effect on sputtering induced 
intermixing and on the observed asymmetry. Hence, we 
expect that the energetic particle scattering mechanism is 
weakened heavily when mass-anisotropy vanishes and the 
magnitude of intermixing is also weakened approaching 
the {R^) values of Ti/Pt. Therefore the observed mix- 
ing asymmetry can solely be explained by the hyperther- 
mic light particle backscattering mechanism at the mass- 
anisotropic interface. The stronger simulated mixing in 
Ti/Pt than in Cu/Co could mostly also be accounted for 
by mass effect. The artificially set mass-isotropy in Ti/Pt 
results in the decrease of (i?^) (filled diamonds in Fig. 
2) which drops to the regime of Cu/Co. In Ti/Pt, the 
heavy atoms at the interface block the transport path of 
the film atoms and reverse the flux of energetic Ti atoms. 
In the artiflcially mass-isotropic Pt/Ti (iso- Pt/Ti), how- 
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FIG. 4. The Cleri-Rosato potential energy as a function of 

the interatomic distance (A) in Ti and in Pt averaged for an 
atom with its 12 first neighbors in the hep and fee lattices. 



ever, (i?,^) is similar to that of in iso-Ti/Pt and somewhat 
larger than in Cu/Co (see Fig 2) hence mass effect does 
not solely explain intermixing in these bilayers. 

To understand the physical origin of the higher inter- 
mixing in iso-Pt/Ti and in iso-Ti/Pt than in Cu/Co, 
we explain in general intermixing by a more complex 
anisotropy of the interface then previously expected. We 
can rule out the effect of thermochemistry (e.g. heat 
of mixing AH) because it has also been shown recently, 
that AH has no apparent effect on intermixing during 
single-ion impact in this bilayer [15]. Also, thermochem- 
istry could not explain the mixing asymmetry observed 
in Ti/Pt and in Pt/Ti since AH must be the same in 
magnitude in the two systems. Moreover, simulated 
ion-sputtering with AH ^ (switching off the attrac- 
tive term of the cross-potential and using purely repul- 
sive (Born-Mayer) Ti-Pt interaction potential for Pt/Ti) 
leads to intermixing instead of the expected sharp inter- 
face (see also animations [14]). We reach the conclusion 
that there must be another anisotropy in the system. 

In order to find another anisotropy paramater, the em- 
ployed interaction potential energy functions are plotted 
in Fig (4) as a function of the interatomic distance for 
Ti and for Pt. It can clearly be seen that the thermal 
atomic size (the width of the potential valley at higher 
atomic kinetic energy) is different. The atomic radius 
(the width of the vibration around the equilibrium lat- 
tice site) of the Ti atoms is larger than that of in Pt. The 
atomic size anisotropy is even more robust at 1-2 eV en- 
ergies above the energy minimum (hyperthermal atomic 
energies, see the horizontal lines in Fig 4). This explains 
the fact that the smaller and heavier Pt atoms are the 
first ballistic diffuser during ion-beam mixing [15]. The 
atomic size difference (with smaller contribution) and the 
mass-anisotropy together could explain the stronger sim- 
ulated (and measured) intermixing in Pt/Ti and in Ti/Pt 
than in Cu/Co. In the latter mass-isotropic bilayer no se- 
rious atomic size difference (anisotropy) is found (see the 
inset Fig (4)). The mixing in the iso-Pt/Ti could also 
be understood by the size anisotropy: the larger ener- 



getic Ti atoms are reversed at the interface and confined 
in the bulk. The smaller energetic Pt atoms, however, 
are injected to the Ti bulk during the cascade period. 
Therefore the BHP mechanism is governed by a complex 
anisotropy of the interface system: both by mass and 
size anisotropy. In the iso-Ti/Pt the energy density of 
the cascade is weaker due to the proximity of the free 
surface hence the influence of the size difference is less 
efficient on intermixing. 

We could explain the observed mixing asymmetry be- 
tween Ti/Pt and Pt/Ti in terms of the backscattering of 
energetic particles at the mass-anisotropic interface. The 
mass effect might not solely explain intermixing in Pt/Ti: 
the atomic size and mass-anisotropy might together be 
responsible for the observed strong interfacial mixing and 
the weaker mixing in Cu/Co is explained by atomic size 
and mass-isotropy. 
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